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Abstract 

The leading paradigm for chaotic flows dominated by dissipation predicts an exponential 
growth of the mean distance between pairs of fluid elements, in the long run. This is reflected 
in the analysis of experimental results on tracer particles and the discussions which follow, as 
reported in recent experimental and numerical publications. To quantitatively validate this 
prediction, we have conducted a nricrofluidic experiment generating elastic turbulence , a flow 
characterised in the literature as smooth in space and random in time. To our great surprise, we 
discovered that the pair separation follows a much slower power-law — also known as ballistic 
a notion overlooked so far for flows of this type. We provide conclusive experimental evidence 
that this scaling is well-desribed by the same coefficients derived from the short-time dynamics. 

Our finding reinforces the role of the ballistic regime over a significant range in time and space, 
providing a quantitative estimation for the spreading of particles in mixing microfluidic flows 
based on the initial velocity snapshot. Finally we note that the conditions for the asymptotic 
exponential pair separation are quite stringent when it comes to tracer particles and are unlikely 
to be realised in wall-bounded flows. Therefore this work raises questions regarding the relevance 
and applicability of the currently leading paradigm. 
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When we stir sugar in a cup of coffee we typically drive the liquid in circles using the tea-spoon, 
yet the flow quickly evolves into a three-dimensional chaotic one, tremendously accelerating the 
homogeneous distribution of the sweetener throughout the beverage. Turbulent flows, typical in 
our everyday life, are renowned for their efficient mixing and intensification of particle dispersion, 
having impact on natural processes as well as wide applications in the industry [1], At their smaller 
scales — the so-called dissipative or Kolomogorov scales — such flows are characterised by a velocity 
field which is smooth in space yet fluctuates strongly in time, also known as the Batchelor regime 
[ 2 ] - 

Understanding transport phenomena at small scales is of importance for the study of combus¬ 
tion, contamination and the deformation of polymers by the flow [3], to name but a few physical 
examples, as well as for the study of biological processes such as cytoplasmic streaming [4], com¬ 
munication via chemo-attractants and food searching [5] . Moreover, it has immediate technological 
implications, mainly when it comes to miniature devices. Microfluidic systems have become a com¬ 
mon tool in research and industry. As a prominent part of lab-on-a-chip apparatuses, they are 
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implemented for micro-chemistry [ 6 , 7], soft condensed matter [ 8 ], bioanalytics such as PCR [9, 10], 
biomedicine and other research and engineering applications [11]. Nevertheless, there is limited 
knowledge of transport phenomena in these devices beyond the scenario of laminar steady flows, 
which is not the only case of interest [ 8 , 12 , 6 ]. 

The study of particle dispersion in flows is at the basis of the understanding of transport 
processes [2, 1, 3]. In the early 1950s, G. K. Batchelor predicted that the mean length of material 
lines in turbulent flows would grow exponentially in the course of time, in the long time limit; this 
stems from the notion that the line elements it consists of can be considered short enough such 
that the distance between the ends of an element remain within the dissipative scale throughout 
the motion [13, 14]. In some recent works the terms material line and material line element are 
used interchangeably to indicate the separation vector between two passive particles in the fluid 
[15, 16, 3]. And indeed Batchelor’s prediction has been later reformulated for tracer particles in the 
form of exponential pair separation; see [3, §2] for example. To illustrate how this comes about, let 
us consider a pair of passive tracers separated by the vector R, and whose relative velocity is u. 
The evolution of the squared separation distance R 2 follows 


1 d_ 

2 dt 


R 2 = u • R = uiR , 


( 1 ) 


where ui is the separation velocity, defined by the above relation. The analysis then proceeds by 
assuming a linear flow approximation — the pair separation is considered to be small enough such 
that the velocity of one tracer is linearly related to that of the other. Using this approximation 
Eq. 1 reduces to 


1 d 

2 dt 


R 2 


£R 2 , 


(2) 


where £ no longer depends on R. For the case of chaotic flows, £ can be modelled as a random 
variable, and analyses often focus on the expectation value of such equations [14, 2, 3], Additionally 
assuming the correlation time of £ to be very short compared to the observation time [2, 3], a 
generalisation of the central limit theorem - the multiplicative ergodic theorem of Oseledec [ 2 , 17] 
- is applied, resulting in the exponential pair dispersion 


{R(t) 2 } = R 2 (t 0 ) exp |2£(f-f 0 )j . 


( 3 ) 


This is a relation for the time evolution of the second moment of pair separation distances. In this 
form one can identify 2£ with the generalised Lyapunov exponent of the second order, which is 
generically not trivially related to the ordinary (maximal) Lyapunov exponent; see [18, §3.2.1], [17, 
§5.3], [2] and others. 

Much of the theoretical and numerical literature discussing pair dispersion in the dissipative 
sub-range is devoted to the evaluation of £ in terms of the typical time-scale of the flow r, that is 
£ = 7 /V; see [13, Eq. 2.9]. The value of 7 is still under debate as can be learnt from [14, §24.5] as 
well as [3] and references therein. 

At the time of writing, the exponential pair dispersion is regarded as the leading paradigm for 
chaotic flows which are spatially smooth, as manifested by the analysis of recent experimental results 
and the discussions which follow. Jullien [19] studied an instance of the Batchelor regime flow in 2d 
turbulence, where the velocity field was inferred experimentally followed by numerical integration 
of tracers simulated on a computer; an exponential separation, referred to in that context as Lin’s 
law [ 20 ], was reported during an intermediate time interval of between one to twice the value of the 
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estimated r, after which a power-law scaling has been observed. Salazar and Collins [3] estimated 
7 from measurements of (£) in 3d turbulence reported by Guala et al. [16]; this estimate should 
be taken with grain of salt not only because it is unclear whether these measurements are indeed 
restricted to the dissipative scales but also as, although related, the quantities (£) and £ are not 
the same [14]. Even more recently, Ni and Xia [21] reported measurements in 3d turbulent thermal 
convection and inferred 7 from exponential fits to the mean squared pair separation distance; as 
presented in [21, Fig. 1], the fits are taken at time intervals of up to one Kolmogorov time-scale, 
a time too short with respect to the underlying assumptions, and thereafter the data grows faster 
than the evaluated exponentials. We therefore find it safe to conclude that the literature on the 
subject is lacking conclusive experimental evidence. 

In this paper we report experimental results of particle pair dispersion where, in contrast to the 
above mentioned publications [19, 3, 21], the memory of the initial relative velocity prevails the 
dynamics, showing no signature of the asymptotic exponential growth. 

Studying pair separation in the dissipative sub-range over long times in intense turbulence poses 
a technological challenge. The high velocities, typical of high Reynolds flows, restrict the length 
of the obtained trajectories as exemplified by the above mentioned reports [19, 15, 21] and other 
recent works [22], This is one of the reasons for which the experimental literature on pair dispersion 
in smooth chaotic flows is lagging behind the theoretical one. 

However, chaotic flows do not necessarily require fast velocities or large vessels. One might 
expect viscous flows to be regular in microfluidic tubes and the question of mixing in these appara¬ 
tuses may seem irrelevant due to their low Reynolds nature. Nevertheless, when a minute amount 
of long flexible polymers, such as DNA and protein filaments, are introduced, the flow may develop 
a series of elastic instabilities which render it irregular and twisted, a phenomenon termed elastic 
turbulence [23, 8 , 24], 

The study of pair separation dynamics in elastic turbulence, taking place inside a tiny tube, 
presents technical challenges: (i) The positions of tracers are needed to be resolved over long 
times and distances, in particular when the tracers get nearby to each other, whereas the flow 
is chaotic and three-dimensional; (ii) The scales at which the dynamics takes place require the 
use of a microscope, where 3d imaging is non-trivial; (iii) The flow fluctuations in time dictate 
a high temporal resolution; (iv) The statistical nature of the problem demands a large sample of 
trajectories, which in turn requires long acquisition times and reliable automation. To satisfy these 
requirements a novel method was implemented. The 3d positions of the fluorescent particles are 
determined from a single camera 2 d imaging, by measuring the diffraction rings generated by the 
out-of-focus particle; this way the particle localisation problem turns into a ring detection problem, 
which is addressed accurately and efficiently in Ref. [25]. 

Below we report the results of an experiment designed to measure the value of 7 using direct 
Lagrangian particle tracking in an elastic turbulence flow in a microfluidic tube. 


Results 

We have established an experimental database of trajectories derived from tracers in the flow; see 
Materials and Methods. One example of a pair realisation can be found in Fig. 1 A,B. A sub-sample 
of pair separation R (St) curves is plotted in Fig. 1C; here St = t — to and Ro = R (St = 0). 
These curves, together with the rest of the pairs which were found in each Ro bin, were collected 
and the second moments were calculated to construct the datasets of ([R(St) / Rq) )r 0 presented 
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in Fig. 2 on a semi-logarithmic scale, where an exponential curve would have appeared linear (the 
un-normalised data (R 2 (5t))R 0 can be found in Fig. S7; the sample size data is presented in Fig. S8). 



Figure 1: Pair dynamics example. The trajectories of two tracers are plotted in the left panels. The right panel 
shows a sub-sample of pair separation distances in the course of time. The figure demonstrates the chaotic nature 
of the flow as manifested by pairs, and outlines the analysis forming the ensembles of pairs. (A) A projection on to 
the camera plane, overlaid on a bright field image of the observation window (further technical details are provided 
in Fig. SI and in the Materials and Methods). ( B ) A side projection. For the sake of visualisation the vertical axis 
is stretched by 3/2. All pairs of tracers which were detected at some instant at a prescribed separation distance, 
Ro = 10 ± 0.5 pm in this example, are collected to form one ensemble. The event at which the pair separation was 
nearest to Ro, marked by the red circles in the plot, is recorded as to for the particular pair for later analysis. The 
colour code in the plot denotes time, which spans 4 s in this case. (C) A sub-sample of pair separation distance 
trajectories R(5t) for pairs belonging to the Ro = 10 pm ensemble, presented on a semi-logarithmic scale; For each 
pair, St = t — to is the time elapsed since to- The colour code denotes time, scaled separately for each curve. The 
plot shows a sub-sample of 49 pairs out of more than 4 x 10 4 pairs in this ensemble; see Fig. S8 for sample size data. 

However, our data shows no supporting evidence for the predicted exponential growth of the 
second moment of pair separations. The insets present a zoom-in on the temporal sub-intervals 
where the full range plot may seem to contain linear segments. Nevertheless, no single slope, to 
be identified with 2£ in Eq. 3, can be found. Moreover, an exponential pair dispersion should 
extrapolate to the origin on this plot, and this is clearly not the case. 

To our great surprise we have found that, for a significant fraction of the observation time, 
the mean relative pair dispersion evolves quadratically in time to leading order (||R(<5f) — Ro|| 2 ) oc 
5t 2 ; this observation is evident in the insets of Fig. 3, and stands in sharp contradiction to the 
interpretations of recent experimental results [19, 3, 21]. To better understand the source for this 
scaling let us write the Taylor expansion around 5t = 0 

R (5t) = Rq + uq 5t + -uq 5t 2 + O (5t 3 ) . (4) 
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Figure 2: Pair dispersion normalised by the initial separation. The plot shows the average squared pair 
separation distance, normalised by the initial separation, {(R(5t) / Ro) )r 0 for various Ro between 6 and 50pm; The 
insets show a zoom-in on the initial and intermediate time intervals. The data shows no supporting evidence for the 
exponential time dependence which follows from Eq. 3. 


Substituting this in the expression for the relative pair dispersion and considering the ensemble 
average over pairs of the same initial separation 

(||R(5t) - Ro|| 2 )r 0 = {u 2 )R 0 ,t 0 ^t 2 + (u ' u )R 0 ,t 0 St 3 + 0(6t 4 ) , (5) 

we find that the leading order term at short times is indeed quadratic in St — the so-called ballistic 
regime. 

This is a universal property which does not require any assumptions on the character of the flow. 
A sign of this scaling beneath the dissipative scale has been observed in simulations of isotropic 
turbulence [26, Fig. 5]. Yet neither quantitative analysis of the coefficients in Eq. 5 nor any contrast 
with the exponential separation prediction were made. Other experimental [27] and numerical 
[26, 28] results are limited to the inertial subrange of turbulence. 

To test this further we rescale the relative pair dispersion by the pre-factor, the mean initial 
squared relative velocity (it 2 )_R 0 ,t 0 - Unlike the case of inertial turbulence, for elastic turbulence there 
are no exact results nor scaling arguments to derive the coefficients appearing in Eq. 5. Therefore 
we extract them from the experimental data; the (u 2 ) R (h t 0 profile as function of Ro is presented 
in Fig. S2A. Indeed we find that our data admits a scaling collapse with no fitting parameters, 
providing a convincing experimental evidence that these observations are well-described by the 
short time expansion of the relative pair dispersion. The data shown in Fig. 3 seem to deviate 
significantly from the quadratic scaling only after 2-3 seconds. 
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Figure 3: Relative pair dispersion. 

(d) Forward in time (||R(<5f) — Ro|| 2 )h 0 for various initial separations (inset) between 6 and 50pm, collapse initially 
on a single curve which follows a power-law St 2 , once rescaled by the average squared relative velocity at the initial 
time, (u 2 )n 0 ,t 0 - A significant deviation from St 2 is noticed after 2-3s, indicating the time beyond which higher order 
terms should be considered. ( B ) Backwards in time relative pair dispersion (||R(— St) — Ro|| 2 }_h 0 for the same initial 
separations (inset), showing the same initial scaling collapse as the forward in time. 
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In order to expose the sub-leading contributions to the initial relative pair dispersion, we sub¬ 
tract the backwards-in-time dynamics from the forward one. This way the time-symmetric terms, 
even powers of St, are eliminated. The result, the time asymmetric contributions presented in Fig. 4, 
shows that indeed initially the next-to-leading order correction follows df 3 and that the curves do 
collapse onto one when rescaled by (u • u), the appropriate coefficient in Eq. 5. The values of 
(u ■ u )_R 0i t 0 were, once again, extracted from the experimental data; see Fig. S4 for the Ro profile. 
However, the deviations from this scaling are noticeable earlier than half a second, much earlier 
than the deviations from the ballistic behaviour discussed above. This hints that the deviation 
observed in Fig. 3 is in fact due to higher order terms, potentially an indication of a transition to 
another regime. The fact that the transition takes place at an earlier time for the larger initial 
separations indicates the potential effects of the vessel size and its geometry, as well as the failure 
of the linear flow approximation; see Fig. S2. 
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Figure 4: Time asymmetric contribution to the relative pair dispersion. Taking the difference between the 
datasets presented in the insets of Fig. 3, |(||R(<5f) — Ro|| 2 — ||R(—<5f) — Ro|| 2 )« 0 exposes the contribution of the 
time-asymmetric elements (odd powers in St), shown in the inset (after sign inversion). These data collapse on a 
single curve which follows 8t 3 when rescaled by (u ■ u).R 0l t 0 . The datasets of Ro < 10pm (grey in the legend) are 
omitted from the main figure due to the scatter of the estimator (see Fig. S4). The plot shows a deviation from St 3 
scaling at times shorter than 400 ms, indicating the dominance of higher order (odd) terms at early times and that 
the 5t 3 term alone does not trivially explain the deviation from St 2 , observed in Fig. 3. 
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Discussion 


In the light of Eq. 5, the ballistic evolution at short times should not be too surprising. And yet, 
to our knowledge it has not been discussed experimentally, nor confronted with the exponential 
separation prediction. Moreover, the ‘memory’ of the initial conditions in our data is longer than 
one may have expected. On the one hand, our observations are consistent with the time-scale 5t* 
derived from the ratio of the coefficients in Eq. 5 (see Fig. S5), as deviations from the 5t 2 scaling 
are expected to be noticeable for 5t < 0.1 St*. This is also manifested in the rescaling of the data by 
5t*, presented in Fig. S6, and to be compared with numerical simulations of inertial turbulence [28, 
Fig. 1], On the other hand, during the time the statistics are dominated by the initial conditions, 
single realisations are not monotonic nor predictable (see Fig. 1), that is to say, our observations 
are not a trivial result of the mean flow alone, which in itself is inhomogeneous. Moreover, during 
one second tracers typically travel distances larger than the width of the tube. Given the stochastic 
nature of the dynamics, the spatial structure of the flow and the limited range of the linear flow 
approximation (see insets in Fig. S2), the Taylor expansion in time Eq. 5 was expected to have a 
rather restrictive radius of convergence. Further analysis is required to clarify whether higher order 
terms are of importance and to estimate the Lagrangian correlation time, which accounts for both 
the temporal and the spatial structure of the flow. Hence, in-depth discussion of the transition 
away from the ballistic regime in our system is left for future work. 

Nonetheless, the time for which the ballistic approximation holds is longer than a second, a 
memory of the initial relative velocities which allows a wide prediction horizon. We are unaware 
of any previous discussion of pair dispersion short time statistics neither in the context of elastic 
turbulence nor for the broader class of wall-bounded viscous mixing flows. The implications are far- 
reaching — knowing the average velocity difference at a given scale, which is an Eulerian property, 
one can determine the degree of dispersion, a Lagrangian property, as was recently discussed in 
[29]. We have demonstrated its predictive power over a significant time interval, which provides a 
new paradigm for chaotic flows dominated by dissipation. 

The exponential growth in Eq. 3 relies on two underlying assumptions [13, 2, 3]: (i) the velocity 
held admits a linear approximation in space throughout the observation time; and (ii) the obser¬ 
vation time is much longer than the correlation time of the velocity gradients. These requirements 
are quite stringent and are not fulfilled by our experiment. To the best of our knowledge these 
assumptions have not been met experimentally for tracer particles so far, and yet the exponential 
growth prediction seems to be the leading paradigm in interpreting experimental results [19, 3, 21]. 
Our estimations for the experimental system presented here indicate these may be possibly relevant 
for Rq < 0.1pm and 5t > 10 s. Examining the experimental parameters reported in recent works 
[1, 27, 22] we find that the trade-off between short correlation times and a large enough dissipation 
scale renders them difficult to reach in intense inertial turbulence, as Batchelor himself anticipated 
[13, §1]; also see [14, §24.5]. 

The lesson in this case highlights the limitations of asymptotic analyses for the descriptions of 
real-world flows — for practical purposes the asymptotic dynamics are not necessarily the generic 
description which should be considered but rather the short-time behaviour, here dictated by the 
Eulerian averages at the initial time. According to our findings the Lyapunov exponents picture 
is not the appropriate one to describe the dynamics of wall-bounded elastic turbulence at scales 
larger than few percent of the vessel size. This raises questions regarding the mechanism for polymer 
stretching once their end-to-end distance reaches few microns length. 

Following our work one may reach the conclusion that the Batchelor prediction [13] is irrelevant. 


By all means, this is not the case! Batchelor considered the evolution of the total length of a material 
line while this work, like the above mentioned reports [19, 3, 21], focused on the shortest distance 
between pairs of tracers starting with a finite initial separation. Therefore the question of what 
natural or experimental systems admit the assumptions underlying the exponential pair dispersion 
Eq. 3 for tracer particles persists. 


Materials and methods 

Methods Summary The work presented here relies on constructing a database of trajectories 
in an elastic turbulence flow. Elastic turbulence is essentially a low Reynolds number and a high 
Weissenberg number phenomenon. The former means the inertial non-linearity of the flow is over¬ 
damped by the viscous dissipation. The latter estimates how dominant is the non-linear coupling 
of the elastic stresses to the spatial gradients of the velocity field compared with the dissipation of 
these stresses via relaxation. This is the leading consideration in the design of the flow cell. 

The Lagrangian trajectories are inferred from passive tracers seeded in the fluid. In order to 
study the dynamics of pairs, the 3d positions of the tracers are needed to be resolved, even when 
tracers get nearby to each other. The requirement of large sample statistics dictates the long 
duration of the experiment, which lasts over days. The fluctuations due to the chaotic nature of 
the flow set the temporal resolution at milliseconds. This leads to a data generation rate of about 
180 GB/h. Hence both the acquisition and the analysis processes are required to be steady and fully 
automated. The 3d positions of the fluorescent particles are determined using 2d single camera 
imaging, by measuring the diffraction rings generated by the out-of-focus particle. This way the 
particle localisation problem turns into a ring detection problem [25]. 

Microfluidic apparatus The experiments were conducted in a microfluidic device, implemented 
in polydimethylsiloxane elastomer by soft lithography, consisting of a curvilinear tube having a 
rectangular cross-section. The depth is measured to be 140 pm, the width is approximately 185 pm 
(see Fig. SI). The geometry consists of a concatenation of 33 co-centric pairs of half circles. The 
working fluid consists of polyacrylamide (MW=1.8 x 10' Da at mass fraction of 80 parts per million) 
in aqueous sugar syrup (1:2 sucrose to d-sorbitol ratio; mass fraction of 78%), seeded with fluorescent 
particles (1 micron Fluoresbrite YG Carboxylate particles, PolySciences Inc.) at number density 
of about 50 tracers in the observation volume. The viscosity of the Newtonian solvent, without 
the polymers, is estimated to be 1100 times larger than water viscosity at 22 °C. This leads to a 
polymer longest relaxation time of r p ~ 100 s [30], which is the longest time scale characterising 
the relaxation of elastic stresses in the solution. The flow was driven by gravity. The maximal 
time-averaged velocity was measured to be roughly u max ~ 250 pm/s. This results in a Reynolds 
number Re < 10 -4 and a global Weissenberg number Wi = t p ~ 250. 

Imaging system The imaging system consists of an inverted fluorescence microscope (IMT-2, 
Olympus), mounted with a Plan-Apochromat 20x/0.8NA objective (Carl Zeiss) and a fluorescence 
filter cube; a Royal-Blue LED (Luxeonstar) served for the fluorophore excitation. A CCD (GX1920, 
Allied Vision Technologies) was mounted via zoom and O.lx c-mount adapters (Vario-Orthomate 
543513 and 543431, Leitz), sampling at 70Hz, 968px x 728px, covering 810pm x 610pm laterally 
and the full depth of the tube. The camera control was based on a modification of the Motmot 
Python camera interface package [31], expanded with a home-made plug-in, to allow real-time 
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image analysis in the RAM [25], recording only the time-lapse positions of the tracers to the hard 
drive. 

Lagrangian particle tracking To construct trajectories, the particle localisation procedure, 
introduced in [25], has to be complemented by a linking algorithm. Here we implemented a kine¬ 
matic model, in which future positions are inferred from the already linked past positions. We 
used the code accompanying [32] as a starting point. The algorithm was rewritten in Python [33], 
generalised to n-dimensions, the kinematic model modified to account for accelerations as well, a 
memory feature was added to account for the occasional loss of tracers, and it was optimised for 
better performance. The procedure accounts for the physical process of particles advected by a 
smooth chaotic flow and for the uncertainties. These arise from the chaotic in time nature of the 
flow (“physical noise”) as well as from localisation and past linking errors (“experimental noise”). 
Finally, natural smoothing cubic splines are applied to smooth-out the experimental noise and es¬ 
timate the velocities and accelerations [34, 35]. The smoothing parameter is chosen automatically, 
where Vapnik’s measure plays the role of the usual generalised cross-validation, adapted from the 
Octave splines package [36]. 

Pairs analysis Within the trajectories database we have identified pairs of tracers which were 
found at some instant at a separation distance close to a prescribed initial separation Rq = 
6 , 7,..., 50 pm, to within 5Rq = ±0.5 pm. The initial time to for a trajectory was set by the 
instant at which the separation distance was closest to Rq. This way, each pair separation trajec¬ 
tory R(5t ) can contribute to an Ro pairs ensemble at most once. See Fig. 1. The number of pairs 
considered in each Ro ensemble is plotted in Fig. S8 as function of 5t. 

Examining the ensemble averages of the relative separation velocity at the initial time (ui}n 0) t 0 , 
presented in Fig. S3, we do not find an indication that our sampling method introduces a bias for 
converging or diverging trajectories, at least for Ro < 22 pm. Such a bias could manifest itself in 
non-vanishing relative separation velocity at these scales. We have yet to resolve the reason for the 
negative values observed at the larger scales. 

Our data supports the linear flow approximation assumption at small enough scales, as indicated 
by the ensemble averages of the initial relative separation velocity; see the inset of Fig. S2 B where 
{uf'jTio.toIRl levels-off at Ro < 12 pm. The same regime is not reached for the relative velocity, yet 
the data in the inset of Fig. S2A does not rule out this possibility for smaller scales. 

One note of warning should made as to the potential limitations when inferring the Eulerian 
mean squared relative velocity from direct Lagrangian particle trajectories, e.g. using (u 2 )r o j 0 
in Fig. S2 as an estimator for the Eulerian second order structure function. The tracer number 
density in our experiment sets the mean separation distance to a value which is larger than the 
initial separations considered in this work. The analysis of the possible consequences are beyond 
the scope of this report. Yet we have no evidence for a discrepancy in the quantities we have 
considered here. 
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Supporting figures 



Figure SI: Microfluidic tube geometry. A bright field image of the observation window; the tube geometry is 
based on a planar concatenation of co-centric pairs of half circles; the observation window encloses two out of 33 in 
total; the arrows, and the yellow broken line indicate the tube walls and their dimensions, specified in microns; the 
tube has a rectangular cross-section of 140pm depth (perpendicular to the image plane). 
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Figure S2: Mean squared relative velocity at to. The averages are taken at to, when the pairs separation 
distance is closest to Ro, as function of Ro, for: (A) squared relative velocity (u 2 )r 0 ^ 0 , which is our estimator for 
the coefficient of the quadratic term in Eq. 5 (see the closing note of the Materials and Methods), and ( B ) squared 
relative separation velocity (M 2 )j? 0l t 0 , where ui = u ■ R/f?. 

The error bars indicate the margin of error based on the 95% confidence interval. The insets show the same data 
rescaled by the squared initial separation distance Ro, plotted vs. Rq. The (w 2 )|i? 0 ,t 0 /Ro data (inset of B) levels off 
as Ro approaches the smaller distances, providing supporting evidence for the linearity of ui with R at scales smaller 
than 12 pm. 
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Figure S3: Mean relative separation velocity at to. The averages are taken at to, when the pairs separation 
distance is closest to Ro, as function of Ro, for the relative separation velocity ( ui)R 0: t 0 ■ The error bars indicate the 
margin of error based on the 95% confidence interval. The inset shows a zoom-in on the smaller range of Ro studied 
here, where (ui)i? 0 ,i 0 is nearly zero. This is an indication of a sampling with no preferential for pair trajectories which 
diverge or converge initially. 
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Figure S4: The coefficient of the cubic term in Eq. 5, (u- u) /(„,*„. The averages are taken at to, when the pairs 
separation distance is closest to Ro, as function of Ro, for (u ■ u }R 0 ,t 0 . It is also an estimator for the time derivative 
of {\u 2 )R 0 ,t 0 - The error bars indicate the margin of error based on the 95% confidence. Due to the extent of the 
margin of error for Ro < 10 we did not use these values for the re-scaling in Fig. 4. 
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Figure S5: Transition time from the ballistic regime. 

= | ( u )Rq-. t 0 j (u ■ u}_R 0 ! t 0 j is plotted as function of the initial separation distance Ro. Deviations of the relative 
dispersions (||R — Ro|| 2 }_b 0 from the St 2 scaling are generically expected to be significant for St > 0.1 St*. 



Figure S6: Dimensionless relative pair dispersion. The short-time statistics of the relative pair dispersion Eq. 5 
is reformulated in a dimensionless form when divided through by (n 2 ) R 0 ,t 0 {St Rg ) 2 . The rescaled data, for various 

initial separations between 6 and 50 pm, is plotted in the Egure against the dimensionless time St j St R(J ; here St Ro 

is the same as plotted in Fig. S5. The datasets indeed collapse onto a single curve {St/St Ro ) 2 — {St/St Ro ) 3 , shown 
by the dashed black line, for St/St Ro < 0.2. As can be seen from the zoom-in (inset), the datasets of larger initial 
separations Ro > 23 pm follow the short-time statisitcs more closely longer. The uncertainty of our estimator for St* Ro 
at smaller intial separations is larger, as reflected by the errorbars in Fig. S5; reducing the degree of the uncertainty 
may provide better understanding of why the data collapse is less pronounced for the smaller Ro ensembles. 
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Figure S7: Pair dispersion. The plot shows the average squared pair separation distance, (R 2 )r 0 for various Ro 
between 6 and 30 pm; the insets show a zoom-in on the initial and intermediate time intervals. This is the same data 
plotted in Fig. 2, only not normalised by Ro- 
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Figure S8: Sample sizes. The number of pairs in the datasets presented in this work plotted against St for the 
various initial separation distances Ro- These sample sizes result from several filters on the observed pairs: each 
single trajectory is required to span more than 140 pm; every R(t) trajectory is let to contribute to an Ro bin at most 
once, at the point where it is nearest to the centre of the bin, this moment is denoted to for the pair in this bin; the 
pair trajectory is required to contain at least 15 measurement points before and after to. Averages of sample sizes 
smaller than 100 were excluded. 
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